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SUMMARY 

The significance of hub flexibility In the nonlinear static and dynamic 
analyses of Advanced Turboprop blades Is assessed. The blade chosen for this 
study Is the 0.175 scale model of the GE-A7-B4 Unducted Fan blade. A procedure 
for coupling the effective hub stiffness matrix to an MSC/NASTRAN finite ele- 
ment model Is defined and verified. A series of nonlinear static and dynamic 
analyses are conducted on the blade for both rigid and flexible hub configura- 
tions. Results Indicate that hub flexibility Is significant In the nonlinear 
static and dynamic analyses of the GE-A7-B4. In order to Insure accuracy In 
analyses of other blades, hub flexibility should always be considered. 


INTRODUCTION 

Because of the potential for very high propulsive efficiency at cruise 
speeds up to Mach 0.8, advanced forms of the propeller, called propfans, are 
being seriously considered for aircraft propulsion. To obtain maximum aero- 
dynamic and acoustic performance, the trend In advanced high speed propeller 
design has been toward thin, swept blades of complicated structural design 
(fig. 1). A research program to establish the required technology for success- 
ful design of propfans Is In progress at the NASA Lewis Research Center (refs. 

1 and 2). 

Part of this effort Is to understand and predict the structural and 
dynamic behavior of these blades. Normally the MSC or COSMIC NASTRAN finite 
element computer programs are used to calculate In-vacuum steady state dis- 
placements (refs. 3 and 4). 

To date, most analyses have assumed that the blades are rigidly attached 
at their base (fig. 2). This simplification makes blade modeling easier, but 
does not allow for the effect of hub flexibility to be Included In the anal- 
ysis. The purpose of this study Is to determine the effect of hub flexibility 
on propfan steady state displacement, frequencies, and mode shapes. 

The effect of hub flexibility Is evaluated by performing steady state 
displacement and frequency analyses on the GE-A7-B4 Unducted Fan (UDF) blade. 
Selection of this blade Is based on the availability of an effective hub stiff- 
ness; characterized by a 6 by 6 stiffness matrix corresponding to the 6 degrees 
of freedom of the hub attachment point. Unfortunately, a finite element model 
of the hub Is not available. To Incorporate this effective hub stiffness Into 
the NASTRAN model, a procedure Is developed for coupling a symmetric stiffness 
matrix to a NASTRAN finite element model. Verification of this procedure, and 
application to the GE-A7-B4 UDF blade. Is presented. 


PROCEDURE 


This section discusses the coupling procedure for connecting the propfan 
blade to the effective hub stiffness matrix using the capabilities on NASTRAN. 
Once this procedure has been Implemented, the coupled blade/hub NASTRAN model 
can be used for both the large displacement (solution 64) and the normal modes 
(solution 63) analyses (ref. 5). The coupling procedure entails, (1) defining 
a grid point at the blade to hub attachment location (fig. 3), (2) connecting 
the grid point defined In 1 to the blade finite element model, and (3) connect- 
ing the grid point to the effective hub stiffness matrix. 

Either an existing grid point or a newly created grid point may be used 
for defining the hub attachment location. If an existing grid point coincides 
with the attachment location, then a new grid point Is not required. Other- 
wise, a new grid point must be created. In figure 3, a new grid point was 
created because none of the existing grid points coincided with the hub attach- 
ment location. 

After the grid point Is defined. It Is attached to Its adjacent grid 
points by rigid bar elements (see ref. 5 for a description of RBAR elements). 

In the figure, four rigid bar elements are shown which connect the new grid 
point to four existing grid points on the blade's finite element model. The 
effect of using the rigid elements Is to constrain the new, and four adjacent, 
grid points to displace as a rigid body. This causes a small region of the 
finite element mesh to be overly stiff. Since only a small region Is con- 
strained, the overall characteristics of the blade will not be affected. 

The effective hub stiffness matrix Is Incorporated Into the blade model 
by using NASTRAN spring elements ( CELAS2 cards). Normally CELAS2 cards 
(fig. 4) are used to connect the degrees of freedom of two unique grid points. 
For this application, the spring elements are used to specify the coefficients 
of the effective hub stiffness matrix. This utilization Involves connecting 
springs from the attachment grid point to ground, and using spring elements to 
specify the coupling between the degrees of freedom at the attachment grid 
point. The CELAS2 spring elements are used for this application because the 
NASTRAN large displacement solution sequence (solution 64), unlike some of the 
other solution sequences, does not permit the direct Input of a stiffness 
matrix. 

Twenty one CELAS2 cards are needed to specify the 36 elements of the 6 by 
6 effective hub stiffness matrix. Six CELAS2 cards are used for connecting 
the attachment grid point to ground In the direction of each of the six global 
degrees of freedom (u x ,Uy,u z ,e x ,e y ,e z ) . The rest of the CELAS2 cards (15 of 
them) are used for specifying the cross coupling between the degrees of freedom 
at the attachment grid point. 

The stiffness coefficients for the six grounded terms are computed by sum- 
ming all of the coefficients In each row of the effective hub stiffness matrix. 
Each stiffness coefficient Is defined by; 

k 1 = K hub 1j = ] » 6 ) 
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The fifteen springs used to couple the degrees of freedom at the attachment 
grid point are defined by; 


k 1j 3 ~ K hub 1j 


where only the off-diagonal coefficients above the diagonal of the hub stiff- 
ness matrix are used. The coupling between the degrees of freedom at the 
attachment grid point Is specified by Inputlng each of the k^j th stiffnesses 
Into field three of the CELAS2 cards. The location of the coupling Is speci- 
fied by Inputlng the I** 1 and degree of freedom numbers Into fields five 

and seven, respectively. The attachment grid point Identification Is specified 
In both fields four and six of the CELAS2 card. 

As an example, assume that It Is desired to connect the stiffness matrix 
shown In figure 5(a) to an attachment grid point "Gl". (A 2 by 2 stiffness 
matrix Is shown for simplicity. In an actual application the stiffness matrix 
will probably be 6 by 6.) The first row (1 =1) of the stiffness matrix Is 
summed (300+ (-100)) and the resulting spring value Is specified In field three 
of the CELAS2 card (fig. 5(b)). A "1" Is specified In field five, and the 
attachment grid point Identification Is specified In field four. Fields six 
and seven are left blank Indicating that the spring Is connected to ground. 

Next, the second row (1 = 2) Is summed and the result Is specified In the 
second CELAS2 card (fig. 5(c)). This spring Is also connected to ground. 

Only two grounded springs are required for this example because the stiffness 
matrix Is only 2 by 2. As previously mentioned, a 6 by 6 matrix requires the 
specification of 6 grounded springs. 

The coupling spring (k^j = -(-100)) Is specified In the CELAS2 card shown 
In figure 5(d). This spring Is used to couple the first and second degrees of 
freedom. The NASTRAN program will automatically Insert a 100 Into the first 
and second diagonal elements of the matrix, and a -100 Into the off-diagonal 
locations. Note that the 100 Inserted on the diagonal Is added to the values 
that were previously specified via the springs In figures 5(b) and (c). For 
this example, only one coupling spring Is required; where as for a 6 by 6 
matrix, 15 coupling springs would be needed. 


VERIFICATION 

The proceeding procedure for connecting the blade to the hub was verified 
with a flat plate consisting of seven brick elements (fig. 6). The model Is 
cantilevered, has a length of 14, a width of 2, and depth of 0.4. The brick 
elements are limited to three translational degrees of freedom at each grid 
point. All seven brick elements have equal size, and material properties. 

Verification of the defined procedure consists of the following: 

(1) A dynamic analysis of the full flat plate model described above. 
Results of this analysis will act as a control for the following analyses. 
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(2) A dynamic analysis of the reduced flat plate model shown In figure 7. 
This model Is generated by substituting element seven with Its effective mass 
and stiffness matrices. These matrices are attached to grid point 33 and then 
point 33 Is attached to the base with rigid elements. The effective mass 
matrix Is applied to grid point 33 with C0NM1 cards. 

(3) A repeat of the analysis In the above omitting the effective mass 
matrix. 

The effective mass and stiffness matrices are derived from element seven's 
property matrices using static condensation (see OMIT card, ref. 5). 

Table I presents the first bending, torsion, and edgewise bending frequen- 
cies recorded for analyses 1 to 3. For the three analyses conducted, both 
bending frequencies compare favorably with, at most, a 2.5 percent frequency 
shift between analyses 1 and 3. Also, It Is apparent that the significance of 
the effective mass matrix Is negligible with at most a 0.03 percent first edge- 
wise bending frequency shift between analyses 2 and 3. However, added tor- 
sional stiffness exists In the reduced flat plate model, resulting In a 
14.3 percent first torsional frequency shift between analyses 1 and 3. As 
mentioned In the proceeding section, this added torsional stiffness In analyses 
2 and 3 Is due to the constraints applied by the rigid elements. These con- 
straints require that the bottom face of element 6 remains quadrilateral, when 
In fact It warps In the torsional mode. Because of the simplified model used 
for verification of procedure, the effects of added torsional stiffness are 
more dramatic. In propfan analysis, It Is safe to assume that the added tor- 
sional stiffness due to constraints will be negligible because the rigid ele- 
ments only constrain a small portion of the blade's base with respect to the 
size of the model, and the number of elements used (fig. 3). 


APPLICATION 

The blade chosen to assess hub flexibility Is the 0.175 scale model of 
the GE-A7-B4. The finite element model of this blade, along with the effective 
stiffness matrix characterizing hub flexibility, has been provided by the 
General Electric Company (fig. 8). To assess the effect of hub flexibility, a 
series of large displacement and dynamic analyses were performed at various 
rotational speeds using MSC/NASTRAN's solution sequence 64 and 63, 
respectively. 

Figure 9 presents the additional cards Included In the bulk data deck of 
GE-A7-B4's finite element model In order to account for hub flexibility. The 
GRID card defines the massless, 6 degrees of freedom grid point, 757, where 
the hub flexibility matrix Is to be applied (fig. 10). The RBAR cards con- 
strain the degrees of freedom of grid point 734, 735, 742, and 743 to the 
degrees of freedom of grid point 757. And the CELAS2 cards add the stiffness 
Influence coefficients, defined by the effective hub stiffness matrix, directly 
Into the global stiffness matrix of the finite element model with spring con- 
stants defined by the procedure section of this paper. 

Results of the analyses are presented In tables II and III, and figures 11 
and 12. Figure 11 presents the magnitude of leading edge tip displacements 
versus rotational speed, for both rigid and flexible hubs. At low rpm, the 
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flexible hub has little effect on the static displacement of the GE-A7-B4. 
However, at rpm's greater than 4000, the Influence of the flexible hub becomes 
notable. Figure 12 presents the first bending, second bending, first tor- 
sional, and third bending eigenvalues for the GE-A7-B4 at various rpm's. In 
both rigid hub and flexible hub configurations. While the first bending and 
first torsional frequencies seem to be unaffected by the flexible hub, there 
Is an appreciable discrepancy between the flexible hub and the rigid hub con- 
figurations concerning the second bending and third bending frequencies. Also 
shown In figure 12 Is an edgewise mode that Is being picked up In the flexible 
hub configuration. The frequency of this mode lies between the second bending 
and first torsional modes. The first edgewise mode was not seen for the rigid 
hub configuration In the lower frequency range. Figure 13 presents nodal line 
plots for the mode shapes discussed above. 


SUMMARY AND CONCLUSION 

The GE-A7-B4 UDF blade was chosen to assess the effects of hub flexibility 
on propfan blade vibrations. This blade was chosen due to the availability of 
a 6 by 6 nodal stiffness matrix characterizing hub flexibility. A finite ele- 
ment model of the hub was not available. 

In order to assess the effects of hub flexibility, a series of large dis- 
placement and dynamic analyses needed to be performed at various rpm's using 
MSC/NASTRAN ' s solution sequences 64 and 63, respectively. While most MSC/ 
NASTRAN solution sequences allow for the modification of the global stiffness 
matrix with stiffness coefficients, there was no direct method to apply the 
hub stiffness matrix to the blade's finite element model In a solution 64 large 
displacement analysis. A procedure was therefore developed, and verified, 
allowing for the coupling between a symmetric stiffness matrix and a NASTRAN 
finite element model In a solution 64 large displacement analysis. Although 
this procedure was demonstrated using an analytically derived effective hub 
stiffness, It Is equally applicable for an experimentally obtained hub 
stiffness . 

The effect of hub flexibility on steady state displacement and frequencies 
has been shown to be significant for the GE-A7-B4 Unducted Fan blade. A drop 
In the second bending, third bending, and first edgewise modal frequencies may 
have a significant Impact on any aeroelastlc analysis conducted on the blade. 

Overall, In order to Insure accuracy In analyses of other propfan blades, 
hub flexibility should always be considered. If hub flexibility effects are 
not assessed, the accuracy of the results will be questionable. 
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TABLE II. - GE-A7-B4 LEADING EDGE TIP 
DISPLACEMENTS: RIGID HUB VERSUS 

FLEXIBLE HUB (PERCENT SPAN LENGTH) 


rpm 

Rigid hub 

Flexible hub 

3000 

0.1878 

0.1931 

4000 

.2922 

.3053 

5000 

.3923 

.4205 

6000 

.4803 

.5348 

7000 

.5527 

.6496 

8000 

.6104 

.7682 
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TABU III. - GE-A/-B4 MODAL FREQUENCIES (Hz) 


rpm 

Mode 

Rigid hub 

Flexible hub 

0 

IB 

198.26 

188.45 


2B 

577.07 

502.77 


IE 


675.21 


IT 

757.46 

748.47 


3B 

1215.19 

1120.14 

2000 

IB 

206.31 

196.18 


2B 

582.83 

508.03 


IE 


678.95 


IT 

764.61 

754.16 


3B 

1220.90 

1127.76 

3000 

IB 

215.89 

205.35 


2B 

589.69 

514.32 


1 C 

• u 


A? 


IT 

773.77 

761.64 


3B 

1227.93 

1137.22 

4000 

IB 

228.49 

217.40 


2B 

598.72 

522.62 


IE 


689.35 


IT 

786.88 

772.59 


3B 

1237.58 

1150.33 

n 

IB 

243.56 

231.78 


2B 

609.45 

532.51 


IE 


696.55 


IT 


787.22 


3B 

1249.67 

1166.88 

6000 

IB 

260.54 

247.94 


2B 

621 .44 

543.58 


IE 


704.86 


IT 

825.12 

805.55 


3B 

1263.99 

1186.58 

7000 

IB 

278.86 

265.47 


2B 

634.26 

555.48 


IE 


714.21 


IT 

849.78 

827.44 


3B 

1280.25 

1209.08 

8000 

IB 

298.38 

283.86 


2B 

647.80 

567.81 


IE 


724.46 


IT 

877.88 

852.37 


3B 

1298.45 

1233.80 


0 





















FIGURE 2. - PRC 


NEWLY CREATED GRID POINT AT 
BLADE-TO-HUB ATTACHMENT 
POINT —7^ 





INPUT DATA CARD CELAS2 SCALAR SPRING PROPERTY AND CONNECTION 


DESCRIPTION : DEFINES A SCALAR SPRING ELEMENT OF THE STRUCTURAL 

MODEL WITHOUT REFERENCE TO A PROPERTY CARD. 


FORMAT 


1 

2 

3 

A 

5 

6 

7 

8 

9 

10 

CELAS2 

EID 

K 

G1 

Cl 

G2 

C2 

GE 

S 













FIELD CONTENTS 

EID UNIQUE ELEMENT IDENTIFICATION NUMBER (INTEGER >0) 

K THE VALUE OF THE SCALAR SPRING (REAL) 

G1,G2 GEOMETRIC GRID POINT IDENTIFICATION NUMBER (INTEGER >0) 

Cl ,C2 COMPONENT NUMBER (6 >: INTEGER > 0) 

GE DAMPING COEFFICIENT (REAL) 

S STRESS COEFFICIENT (REAL) 

FIGURE A. - CELAS2 CARD. 


1 2 

1 T 300 -100 1 

2 L -100 250 J 

EFFECTIVE STIFFNESS MATRIX (2 DOF) 


j~ 200 j 

CELAS2,"EIDr,200,"Gr,1 


200 — 1 

— 150 J 

CELAS2,"EID2'M50, rt G1",2 


(200 + 100 ) -100 ] 

(d) 

L -100 (150+ 100) J 

CELAS2,"EID3'\100,"G1",1,"G1",2 

FIGURE 5. - CONNECTING A 2 DOF EFFECTIVE STIFF- 
NESS MATRIX TO GRID POINT "G1". 
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FIGURE 6. - FLAT PLATE MODEL. 


FIGURE 7. - REDUCED FLAT PLATE MODEL 













\ \ \~~T u\\ 

tUTTlV V 



GRID 

RBAR 

RBAR 

RBAR 

RBAR 

CELAS2 

CELAS2 

CELAS2 

CELAS2 

CELAS2 

CELAS2 

CELAS2 

CELAS2 

CELAS2 

CELAS2 

CELAS2 

CELAS2 

CELAS2 


757 

100 

200 

300 

400 


757 

757 

757 

757 


11 329209. 

22 866345. 

33 40591 . 

44 223558. 
55 243175. 
66 552796. 

12 -55. 

13 -7672. 

14 -15. 

15-125348. 
16 -294. 

23 280. 
24-166324. 


- 0.0102 

742 

743 
735 
734 
757 
757 
757 
757 
757 
757 
757 
757 
757 
757 
757 
757 
757 


Ik] 


HUB 


0.0578 

123456 

123456 

123456 

123456 

1 

2 

3 

4 

5 

6 
1 
1 
1 
1 
1 
2 
2 


*1.958 E5 

5.500 El 

7.672 E3 1.500 tl 

1.253 t3 

2.940 t?~ 

r*~\ uro 
ULLMJ2. 

or 

G7. 

757 

o 





CELAS2 

26 

47071. 

757 

2 

5.500 El 

7.472 E5 

-2.800 E2 1.663 E5 

6.700 El 

-4.707 E4 

CELAS2 

34 

49. 

757 

3 

7.672 E3 

-2.800 E2 

2.478 E4 -4.900 El 

7.935 E3 

5.330 E2 

CELAS2 

35 

-7935. 

757 

3 

1.500 El 

1.663 E5 

-4.900 El 7.085 E4 

-4.200 E2 

-1.316 E4 

CELAS2 

36 

-533. 

757 

3 

7.935 E3 -4.200 E2 

1.095 E5 

6.S90 E2 

CELAS2 

45 

420. 

757 

4 

1.253 E5 

6.700 El 

CELAS2 

46 

13163. 

757 

4 

2.940 E2 

-4.707 E4 

5.330 E2 -1.316 E4 

6.990 E2 

6.115 E5. 

CELAS2 

56 

-699. 

757 

5 


4.9572 


757 

757 

757 

757 

757 

757 

757 

757 

757 

757 

757 

757 

757 

757 

757 


123 

123 

123 

123 


2 

3 

4 

5 

6 

3 

4 

c 

6 

4 

5 

6 

5 

6 
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FIGURE 8. - FINITE ELEMENT MODEL OF GE-A7-B4 AND EFFECTIVE HUB 
STIFFNESS MATRIX. 


FIGURE 9. - MSC/NASTRAN BULK DATA DECK CARDS ACCOUNTING FOR HUB 
FLEXIBILITY. 




FIGURE 11. - GE-A7-B4 LEADING EDGE TIP DISPLACEMENTS. RIGID 
HUB VERSUS FLEXIBLE HUB. 



FIGURE 13. - GE-A7-B4 NODAL LINE PLOTS. 
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